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Abstract: Abnormal expression of micro-ribonucleic acid (miRNA) might be clinically valuable 
as a biomarker or treatment target in the early diagnosis, treatment, and prognosis of tumors. 
However, little is known concerning abnormal miRNA expression of laryngeal carcinoma, one 
of the most commonly encountered head and neck tumors. Microarray analysis was used to 
obtain miRNA-expression profiles of ten pairs of freshly frozen laryngeal carcinoma tissue and 
surrounding normal tissue specimens. Characteristic miRNAs that were significantly related to 
laryngeal carcinoma were identified. Verification was performed using an additional 32 pairs 
of samples. The expression of two miRNAs (miR-21-3p and miR-106b-3p) was upregulated 
in both microarray and quantitative real-time polymerase chain-reaction analyses, whereas the 
expression of six miRNAs (let-7f-5p, miR-10a-5p, miR-125a-5p, miR-144-3p, miR-195-5p, 
and miR-203) was downregulated. The decreased expression of let-7f-5p and miR-195-5p is a 
novel finding in head and neck cancer. The target genes of these miRNAs were also predicted 
through multiple software programs. The differential expression of miRNAs might be related to 
the early onset and development of laryngeal carcinoma, and may be exploited as new biomark- 
ers and therapeutic targets in the treatment of laryngeal carcinoma. 

Keywords: microRNA, laryngeal carcinoma, expression profiling, bioinformatics, target 
prediction 

Introduction 

Head and neck tumors constitute one of the main cancers endangering human health, 
ranking eighth in incidence among all cancer types. Laryngeal carcinoma is one of 
the most common types of head and neck tumors, and is the second-most prevalent 
of the respiratory tract tumors. 1 Most laryngeal carcinoma cases are squamous cell 
carcinoma. The incidence of laryngeal carcinoma has been increasing in recent 
years. In the People's Republic of China, the incidence of laryngeal carcinoma has 
increased by about 25% annually in an era of rapid industrialization. 2 3 Although new 
surgical procedures, chemotherapeutic drugs, and more advanced radiotherapy have 
been applied to the treatment of laryngeal carcinoma in the past 30 years, the overall 
survival of patients with laryngeal carcinoma and other head and neck tumors alike 
has not improved, and in some cases has reduced. 4 The survival rate of late-stage 
laryngeal carcinoma is only 30%-40%. 5 Possible epidemiological^ related factors 
for laryngeal carcinoma include smoking, drinking, air pollution, and occupational 
factors. 6 7 The pathogenesis of laryngeal carcinoma is related to the activation of 
oncogenes, including BCL2 and c-MYC, 8 - 9 and the inactivation of tumor-suppressor 
genes, including p53 and RB. 10 n However, the pathogenic mechanism of laryngeal 
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carcinoma at the molecular level is still unclear, which has 
hindered therapeutic interventions and efforts to increase 
the survival rate. 

Micro-ribonucleic acids (miRNAs) are small 
(19-25 nucleotides) noncoding RNAs. 12 At the post- 
transcriptional level, miRNAs regulate the expression 
level of target genes, which influences the regulation of 
biological events, including individual development, apop- 
tosis, proliferation, and differentiation. 1314 Recent studies 
have implicated miRNAs in the regulation of oncogenes 
or tumor-suppressor genes that are germane in the occur- 
rence, development, invasion, and metastasis of a variety of 
tumors. 1516 Furthermore, miRNAs might be related to the 
pathogenesis of head and neck tumors, such as laryngeal 
carcinoma. 17 The differential expression of miRNAs might 
have potential clinical value as biomarkers or treatment 
targets in the early diagnosis, treatment, and prognosis of 
laryngeal carcinoma. 

Although miRNA-expression profiling of many head and 
neck tumors, including oral and thyroid carcinoma, has been 
reported 18 ' 19 the situation is far less clear for laryngeal carci- 
noma, one of the most common head and neck tumors. The 
present study was undertaken to address this shortcoming. 
The aim of this study was to identify new diagnostic 
biomarkers for the early diagnosis of laryngeal carcinoma or 
potential treatment targets for laryngeal carcinoma. 

Materials and methods 

Tissue sampling and clinical 
data collection 

A total of 42 pairs of surgically removed laryngeal car- 
cinoma tissue and paracancerous tissue samples were 
collected from patients who received surgical resection treat- 
ment at Guangdong General Hospital from October 2007 to 
March 2009. Informed consent for tissue use was provided 
beforehand by all patients, and the study was approved 
by the ethics committee of Guangdong General Hospital. 
Laryngeal carcinoma was confirmed by pathological 
study. Normal laryngeal mucosa specimens were retrieved 
10 mm outside the negative margin. Samples were placed in 
RNAlater TissueProtect Tubes (Qiagen, Hilden, Germany) 
and stored at -80°C. Postsurgical pathology confirmed 
the diagnosis of squamous cell carcinoma in all patients. 
The patients did not receive radio- or chemotherapy before 
surgery. Samples were randomly allocated into two groups 
for analysis using the miRNA microarray (ten pairs) 
or quantitative real-time polymerase chain reaction 
(qRT-PCR; 32 pairs). 



RNA extraction 

Total RNA was isolated using TRIzol® (Life Technologies, 
Carlsbad CA, USA) and miRNeasy mini kit (Qiagen) accord- 
ing to the manufacturers' instructions, which efficiently recov- 
ered all RNA species, including miRNAs. RNA quality and 
quantity was measured using an ND- 1 000 spectrophotometer 
(NanoDrop Technologies, Wilmington, DE, USA), and RNA 
integrity was determined by gel electrophoresis. 

Microarray analysis of miRNA 

Total RNA of ten pairs of laryngeal carcinoma tissue sam- 
ples was randomly selected and microarray analysis was 
completed with the sixth generation of miRCURY LNA™ 
(version 16.0) by Exiqon (Vedbaek, Denmark). 

RNA labeling and array hybridization 

After RNA isolation, the miRCURY Hy3™/Hy5™ power 
labeling kit (Exiqon) was used according to the manufac- 
turer's guidelines for miRNA labeling. After stopping the 
labeling procedure, the Hy3-labeled samples were hybridized 
on the miRCURY LNA array according to the instructions 
in the array manual. Following hybridization, the slides were 
washed several times using a wash buffer kit (Exiqon) and 
dried by centrifugation for 5 minutes at 400 rpm. The slides 
were scanned using a GenePix 4000B microarray scanner 
(Molecular Devices, Sunnyvale, CA, USA). 

Data analyses 

Scanned images were imported into the GenePix Pro 
6.0 software (Molecular Devices) for grid alignment and 
data extraction. Hierarchical clustering was performed 
using MEV (multiexperiment viewer) version 4.6 (TIGR, 
Rockville, MD, USA). 

Validation using qRT-PCR 

Total RNA was extracted from 32 pairs of laryngeal 
carcinoma tissue and paracancerous tissue samples. The 
SYBR green assays (Life Technologies) with U6 as internal 
control was used and qRT-PCR was performed to validate 
the results. 

Reverse transcription 

After total RNA extraction, complementary deoxyribonucleic 
acid was synthesized using 2 ng to 2 ug of total RNA by 
Moloney Murine Leukemia Virus Reverse Transcriptase 
(M-MLV RT) according to the manufacturer's instructions 
(Life Technologies). Stem-loop reverse-transcriptase primer 
was added to total RNA (synthesized by RiboBio, Guangzhou, 
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People's Republic of China). The volume was expanded to 
24 uL with ribonuclease (RNase)-free water and mixed well. 
After incubation at 70°C for 10 minutes, the mixture was 
placed in an ice bath for 2 minutes. The following reagents 
were added to the mixture in sequence: 2 \iL of RNase inhibi- 
tor, 8 uL of 5x M-MLV RT buffer, 2 uL of M-MLV RT reverse 
transcriptase (RNase FF), and 2 \xL of dNTP mixture (Life 
Technologies) and RNase-free water (TaKaRa Biotechnology, 
Dalian, People's Republic of China) to a final volume of 40 uL. 
The mixture was placed at 30°C for 10 minutes, followed by 
42°C incubation for 1 hour, and then the mixture was placed 
at 70°C for 15 minutes. The product of reverse transcription 
was used as the template of PCR. 

qPCR 

Amplification was performed on a fluorescent qPCR instru- 
ment (Bio-Rad Laboratories, Hercules, CA, USA). For each 
reaction, three duplicates were prepared. The total reaction 
volume was 20 |iL, which included 10 uL of 2x Mix SYBR 
Green I (Life Technologies) fluorescent solution, 0.25 uL 
each of forward and reverse primers (10 uM), 1 |iL of sample 
template, and RNase-free water that brought the volume to 
20 uL. The reaction cycle was predenaturation at 95°C for 20 
seconds, followed by 40 cycles of denaturation at 95°C for 10 
seconds, annealing at 60°C for 20 seconds and extension at 
70°C for 10 seconds. An amplification curve was plotted. A 
fluorescence signal was collected during the heating process 
from 60°C to 95°C after 40 cycles, and a melting curve was 
plotted. After each reaction, the threshold cycle (Ct) was cal- 
culated. Relative quantitation of gene expression was analyzed 
with the 2~ AACt method using the obtained Ct values. 

Statistical analyses 

The results of the miRNA microarray were analyzed 
with significance analysis of microarrays (SAM) soft- 
ware. The g-value representing the false discovery rate 
(FDR) was set at 0.05, and miRNAs showing differential 
expression in laryngeal carcinoma tissue and paracan- 
cerous tissue were screened and identified. The results 
of qRT-PCR were processed with the 2~ AACt method, and 
the difference between two groups was analyzed with 
the independent-samples /-test using SPSS 13.0 (SPSS, 
Chicago, IL, USA). 

Results 

Demographic and clinical data 
of the 42 laryngeal cancer patients 

The results are summarized in Table 1 . 



Table I Demographic and clinical features of the laryngeal cancer 



patients 


Characteristic 


n (%) 


Age (years), mean ± standard deviation 


S9.8±l 1.7 


Sex 




Male 


39 (92.9) 


Female 


3(7.1) 


Type 




Supraglottic 


12 (28.6) 


Glottic 


28 (66.7) 


Subglottic 


2 (4.8) 


Clinical stage 




1 


19 (45.2) 


II 


8 (19.0) 


III 


2 (4.8) 


IV 


13 (31.0) 


T stage 




Tl 


19 (45.2) 


T2 


10 (23.8) 


T3 


5(11.9) 


T4 


8 (19.1) 


N stage 




NO 


33 (78.6) 


NI-3 


9 (21.4) 


Pathological differentiation 




Poorly differentiated 


3(7.1) 


Moderately differentiated 


23 (54.8) 


Highly differentiated 


16 (38.1) 


Abbreviations: n, number of patients; T. tumor size; N, lymph node status. 

Results of miRNA microarray 



The miRNA microarray analysis screened out 780 miRNAs. 
After SAM analysis with FDR at 0.05, expressions of eleven 
miRNAs were upregulated in tumor tissue and 114 were 
downregulated, compared to expression in normal mucosa 
tissue (Tables 2, 3 and Figure 1). 



qRT-PCR validation 

Eight miRNAs (miR-21-3p, miR-106b-3p, let-7f-5p, 
miR-10a-5p, miR-125a-5p, miR-144-3p, miR-195-5p, and 



Table 2 Significantly upregulated microRNAs (SAM, FDR <5%) 



Name 


Fold change 


q-value (%) 


hsa-miR-2 1 -3p 


1.823432255 


0 


hsa-miR-205-3p 


1.942721294 


0 


hsa-miR-l06b-3p 


1.213735489 


0 


hsa-let-7f- 1 -3p 


1.0752001 18 


0 


hsa-miR-1246 


185.9579442 


I.I 53846 154 


hsa-miR-49l-3p 


4,835.799572 


I.I 53846 154 


hsa-miR-l30b-3p 


1.172966978 


1.978021978 


hsa-miR-3182 


1.816258072 


2.492307692 


hsa-miR-1290 


157.9054397 


2.492307692 


hsa-miR-2 1 l5-3p 


1.252350132 


2.492307692 


hsa-miR-720 


44,171.33165 


4.984615385 



Abbreviations: RNA, ribonucleic acid; SAM, significance analysis of microarrays; 
FDR, false discovery rate. 



OncoTargets and Therapy 20 1 4:7 



submit your manuscript | www.dovepress.com 
Dovepress 



527 



Lu et al 



Dovepress 



Table 3 Significantly downregulated microRNAs (SAM, FDR <5%) 


Table 3 (Continued) 






Name 


Fold change 


q-value (%) 


Name 


Fold change 


„ ,, -~ | , , r . /0/\ 

cj-vaiue \/o) 


hsa-miR-99a-5p 


0.229416618 


0 


hsa-miR-4321 


0.973550709 


1.978021978 


hsa-miR-l228-5p 


0.969412922 


0 


hsa-miR-335-5p 


0.832271437 


1.978021978 


hsvl-miR-H8 


0.822813522 


0 


hsa-miR-l25b-l-3p 


0.95000675 


1.978021978 


hsa-miR-139-Sp 


0.745290916 


0 


hsa-miR-378c 


0.787146273 


1.978021978 


hsa-miR-99a-3p 


0.943420998 


0 


hsa-miR-1243 


0.860868473 


2.492307692 


hsa-miR-l25b-5p 


0.000504141 


0 


hsa-miR-l9l-5p 


0.624357656 


2.492307692 


hsa-miR-208b 


0.856998652 


0 


hsa-miR-20b-5p 


0.810769473 


2.492307692 


hsa-let-7c 


0.163369241 


0 


hsa-miR-299-3p 


0.910445284 


2.492307692 


hsa-miR-26a-5p 


0.41283036 


0 


hsa-miR-IOb-3p 


0.990773272 


2.760467381 


hsa-miR-IOI-3p 


0.01 1547838 


0 


hsa-miR-l45-5p 


0.42142171 1 


2.760467381 


hsa-miR-26b-Sp 


0.044496706 


0 


hsa-let-7g-5p 


0.1 13667898 


2.760467381 


hsa-let-7a-Sp 


0.1 15018014 


0 


hsa-miR-23b-5p 


0.977354147 


2.760467381 


hsa-miR-30a-5p 


0.555963217 


0 


mcv-miR-MI-3p 


0.895527792 


2.760467381 


hsa-miR-IOa-5p 


0.730429293 


0 


hsa-miR-29b-2-5p 


0.988923772 


2.760467381 


hsa-miR-100-Sp 


0.395246953 


0 


hsa-miR-3176 


0.934463957 


2.760467381 


hsa-miR-1297 


0.675961512 


0 


hsa-miR-299-5p 


0.554176248 


2.760467381 


hsa-miR-378a-Sp 


0.959315913 


0 


hsa-miR-3934 


0.965320293 


2.760467381 


hsa-miR-3656 


0.201641898 


0 


hsa-miR-43 1 7 


0.570519513 


2.760467381 


hsa-miR-337-Sp 


0.949658024 


0 


hsa-miR-l95-5p 


0.434556198 


2.760467381 


hsa-miR-375 


0.536732979 


0 


hsa-miR-629-5p 


0.921854174 


2.760467381 


hsa-miR-338-3p 


0.651622173 


0 


hsa-miR-4796-3p 


0.987307839 


2.760467381 


hsa-miR-45 1 a 


3.05E-05 


0 


hsa-miR-363-3p 


0.953261477 


2.760467381 


hsa-miR-3185 


0.985728913 


0 


hsa-miR-208a 


0.952193195 


2.760467381 


hsa-miR-l26-3p 


0.019184657 


0 


hsvl-miR-H6-5p 


0.955520715 


2.760467381 


hsa-miR-497-Sp 


0.83274609 


0 


hsa-miR-3065-5p 


0.984253829 


2.760467381 


hsa-miR-IOI-Sp 


0.970538127 


0 


hsa-miR-1287 


0.98230182 


2.760467381 


hsa-miR-29c-3p 


0.141944203 


0 


hsa-miR-l43-3p 


0.023686522 


3.279352227 


hsa-miR-126-Sp 


0.332925548 


0 


hsa-miR-3650 


0.953404765 


3.279352227 


hsa-miR-625-Sp 


0.839062834 


0 


hsa-miRPIus-1107* 


0.928817452 


3.279352227 


hsa-miR-IOb-Sp 


0.450653681 


0 


hcmv-miR-ULI48D 


0.884013565 


3.279352227 


hsa-miR-30c-Sp 


0.226168321 


0 


hsa-miR-20a-5p 


0.51987351 1 


3.279352227 


hsa-miR-l35a-Sp 


0.93013101 


0 


hsa-miR-l35b-5p 


0.88299591 1 


3.279352227 


hsa-miR-l25a-Sp 


0.09744404 


0 


hsa-miR-200a-5p 


0.982126485 


3.279352227 


hsa-miRPIus-CI070 


0.732527419 


0 


hsa-miR-200a-3p 


0.125234644 


3.279352227 


hsa-miR-l99b-5p 


0.286949641 


0 


hsa-miR-660-5p 


0.919752055 


3.279352227 


hsa-miR-29a-5p 


0.91 1032396 


0 


hsa-miR-3l26-3p 


0.979864468 


3.279352227 


hsa-miR-323b-5p 


0.772414981 


0 


hsa-miR-648 


0.968124533 


3.279352227 


hsa-miR-29c-Sp 


0.862849161 


1.4 160839 16 


hsa-miR-3620 


0.967574384 


4.984615385 


hsa-miR-608 


0.982843666 


1.4 160839 16 


hsa-miR-l30a-3p 


0.473237407 


4.984615385 


hsa-miR-4328 


0.313675459 


1.4 160839 16 


hsa-miR-604 


0.975573233 


4.984615385 


hsa-miR-l29-l-3p 


0.937204791 


1.4 160839 16 


hsa-miR-l38-2-3p 


0.82691 1224 


4.984615385 


hsa-miR-3133 


0.940025146 


1.4 160839 16 


hsa-miRPIus-1137* 


0.909617132 


4.984615385 


hsa-miR-23c 


0.940121803 


1.4 160839 16 


hsa-miR-2053 


0.984695021 


4.984615385 


hsa-miR-30e-Sp 


0.359399404 


1.4 160839 16 


hsa-miR-374c-5p 


0.901568699 


4.984615385 


hsa-let-7f-Sp 


0.653879295 


1.4 160839 16 


hsa-miR-l48a-5p 


0.967409666 


4.984615385 


hsa-miR-30e-3p 


0.734655027 


1.4 160839 16 


hsa-miR-659-3p 


0.92450892 


4.984615385 


hsa-miR-l44-3p 


0.039732846 


1.4 160839 16 


hsa-miR-27b-3p 


0.176138633 


4.984615385 


hsa-miR-98 


0.728012549 


1.4 160839 16 


hsa-miR-652-3p 


0.944030002 


4.984615385 


hsa-miR-203 


0.001955606 


1.4 160839 16 


hsa-miR-557 


0.944551902 


4.984615385 


hsa-miR-486-Sp 


0.71468291 1 


1.4 160839 16 


hsa-miR-34a-3p 


0.978989503 


4.984615385 


hsa-miR-3607-5p 


0.839520455 


1.4 160839 16 


hsa-miR-7-l-3p 


0.992453904 


4.984615385 


hsa-miR-3195 


0.367848181 


1.978021978 


hsa-miR-376a-5p 


0.97514717 


4.984615385 


nsa-miK-zVa-Jp 


U.UU3-<t7/74 1 


1 OTOfn 1 Q"7Q 
1 .7/oUZ 1 7/0 


hsa-miR-l48a-3p 


0.568972471 


4.984615385 


hsa-miR-378a-3p/hsa- 


0.680647548 


1.978021978 


hsa-miR-28-5p 


0.932157662 


t.70to 1 J J 00 


miR-378c/hsa-miR-378d 






hsa-miR-23b-3p 


0.002949212 


4.984615385 


hsa-miR-374a-Sp 


0.629791599 


1.978021978 


hsa-miR-55lb-3p 


0.976024864 


4.984615385 


hsa-miR-30b-Sp 


0.195820835 


1.978021978 


hsa-miR-559 


0.917545663 


4.984615385 


hsa-miRPIus-CI 100 


0.969513061 


1.978021978 














Abbreviations: RNA, ribonucleic acid; SAM, significance analysis 


of microarrays; 






(Continued) 


FDR, false discovery rate. 
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Figure I Cluster-analysis diagram of laryngeal carcinoma micro-ribonucleic acid (miRNA) genes. The heat map of hierarchical clustering shows an miRNA gene differential 
expression profiling between paired laryngeal tumor and surrounding normal tissue samples. Red represents upregulated miRNAs, and green represents downregulated 
miRNAs. 

Abbreviations: It, laryngeal tumor; Int, laryngeal normal tissue. 



miR-203) were validated with qRT-PCR. The expression 
difference of these eight miRNAs was statistically sig- 
nificant between laryngeal carcinoma and normal mucosa 
tissue in the 32 pairs of specimens (P<0.05). Further 
comparison between the microarray and qRT-PCR results 
(Figure 2) showed that the expression pattern of these 
eight miRNAs in microarray was consistent with that in 
qRT-PCR. 

Prediction of target genes 

Potential target genes of the eight miRNAs were prelimi- 
narily predicted using bioinformatics software, including 
PicTar, miRanda, and TargetScan (Table 4 and Figure 3). 

Discussion 

In this study, the latest version of microarray mercury LNA 
(version 16.0) was used to analyze ten pairs of freshly frozen 
laryngeal carcinoma tissue andparacancerous tissue samples, 
and the miRNA-expression profile of laryngeal carcinoma 
was obtained. The results were verified by qRT-PCR using 
another 32 pairs of samples. The results between the two 



analyses were comparable. miR-21-3p and miR-106b-3p 
were upregulated and let-7f-5p, miR-10a-5p, miR-125a-5p, 
miR-144-3p, miR-195-5p, and miR-203 were downregulated 
in the microarray and qRT-PCR analyses. The target genes of 
these miRNAs were also predicted through multiple software 
programs in this study. 

Studies of miRNA-expression profiling of laryngeal 
carcinoma have suffered from some shortcomings. Most 
studies were conducted on head and neck squamous cell 
carcinoma, 17 ' 20,21 and few studies focused simply on laryngeal 
carcinoma. 22 Although the disease locations of all kinds of 
tumors in the head and neck area are in close proximity, the 
larynx is part of the respiratory system, and so is different 
from such components of digestive system as the oral and 
pharyngeal areas. 23 ' 24 Laryngeal carcinoma has signifi- 
cantly distinctive pathological characteristics and prognosis 
from oral cancer, oropharyngeal cancer, and hypopharyn- 
geal carcinoma. For example, the incidence of laryngeal 
carcinoma in men is significantly higher than that of other 
head and neck tumors. 24 In addition, a comparative genomic 
hybridization study showed that the chromosomal patterns 
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miR-21-3p miR-106b-3p let-7f-5p miR-10a-5p miR-125a-5p miR-144-3p miR-195-5p miR-203 

MicroRNA ID 



Figure 2 Comparison of quantitative real-time polymerase chain-reaction (qRT-PCR) results and microarray results of eight micro-ribonucleic acids. The open bars 
represent qRT-PCR analysis, while the solid bars represent the microarray data. 
Abbreviation: RNA, ribonucleic acid. 



of laryngeal carcinoma and other head and neck squamous 
cell carcinomas were significantly different. 25 Thus, without 
distinguishing laryngeal carcinoma from other head and 



Table 4 Putative target genes of differentially expressed micro- 
ribonucleic acids 



MicroRNA 


Chromosome 


Putative target 




location 




miR-2l-3p 


I7q23.l 


LUM, CHD7, NFAT5, WWPI, TGFBI 


miR-l06b-3p 


7q22. 1 


TNFRSF2I, ACPL2, EPHA5, RAB22A, 






SUTRK3 


let-7f-5p 


9q22.32 


ABCB9, CCNJ, RDHIO, GATM, HIC2 


miR-IOa-Sp 


I7q2l.32 


TNRC6B, BCL6, WNK3, RAP2A, GATA6 


miR-l25a-5p 


1 9q 1 3.4 1 


SMARCD2, BMF, ERBB4, MTFI.ZFYVEl 


miR-l44-3p 


I7ql I.2 


ZNF800, ACBD3, TRIO, ETSI, FBXL3 


miR-l95-5p 


I7pl3.l 


TLE4, ST8SIA3, LRP6, WNT3A, KIF2IA 


miR-203 


I4q32.33 


MAP3K13, ZNF28I, FAT3, KIAA1429, 






ADAMTS6 



Abbreviation: RNA, ribonucleic acid. 



neck squamous cell carcinoma, the results will be unavoid- 
ably biased. 

In this study, the differentially expressed miRNAs 
detected in laryngeal carcinoma were different from reports 
of other head and neck tumors or laryngeal carcinoma. 
The numbers of differentially expressed miRNAs were 
larger, and were predominantly accounted for by down- 
regulated miRNAs. The sample size of laryngeal carcinoma 
samples in this study was the largest among studies on head 
and neck tumors reported to date. Due to the heterogene- 
ity of tumor samples, a large sample size might be ben- 
eficial for the accurate detection of abnormally expressed 
miRNA. 26 In addition, the present study utilized the latest 
version of microarray chip with locked nucleic acid (LNA) 
technology, 27 which guaranteed the maximal degree of reli- 
ability of the results. This is the first report on data from 
southern People's Republic of China, and the different geo- 
graphic locales of patients might also have contributed to 
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Figure 3 Target-gene prediction of selected micro-ribonucleic acids (miRNAs) through bioinformatic techniques. Red-box nodes represent miRNA, and green-cycle nodes 
represent putative target gene. Edges describe the inhibitory effects of miRNA on potential target gene. 



the difference of the results. Furthermore, some studies only 
performed microarray detection without further verification 
by qRT-PCR, which has been recognized as a gold standard. 
This might also affect study results. 

Among the eight differentially expressed miRNAs in 
laryngeal carcinoma verified by qRT-PCR, miR-21-3p, 22 
miR-106b-3p, 20 miR-10a-5p, 20 miR-125a-5p, 17 miR-203, 28 



and miR-144-3p 29 have been reported in other head and neck 
cancer studies. miR-125a-5p, miR-203, and miR-144-3p, 
which were only reported in other studies utilizing microar- 
rays, were further verified by qRT-PCR in our study. 

In addition, miRNAs not previously reported or inconsis- 
tent with previous head and neck squamous cell carcinoma 
or laryngeal carcinoma studies were also identified and 
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included let-7f-5p and miR-195-5p. Lethal-7 (let-7) was 
initially discovered in Caenorhabditis elegans. Being one 
of the first tumor-suppressive miRNAs identified, let-7 is 
downregulated in lung, ovarian, and breast cancer. Studies 
showed that let-7i is highly expressed 17 in head and neck 
squamous cell carcinoma, and the expressions of let-7a, c, 
and e are low. 20 In our study, let-7f-5p showed low expres- 
sion, which was opposite to a report from Tran et al. 30 The 
samples used in the prior study were from the cell lines of 
different head and neck squamous cell carcinomas, includ- 
ing tongue, tonsil, hypopharynx, and larynx, and no further 
verification of let-7f was conducted. 30 Our study used freshly 
frozen laryngeal carcinoma tissue as samples, and qRT- 
PCR verification was performed. This further confirmed 
the advantage of singling out such diseases as laryngeal 
carcinoma, instead of studying a mixture of different types 
of head and neck tumors. miR-195 belongs to the miR-15/ 
miR-16 family. miR-195 expression was downregulated 
in some tumor types, including hepatocellular carcinoma 
and breast cancer, 31 suggesting that MIR195 might be an 
important tumor-suppressor gene. Also, it has been reported 
to be involved in the cell cycle, apoptosis, and proliferation 
through regulating such target genes as WEE1, CDK6, and 
BCL2. M Our study is the first to demonstrate the possible 
presence of abnormally decreased expression of miR-195 
in laryngeal carcinoma. 

The molecular mechanism of miRNA in the regula- 
tion of laryngeal carcinoma-related genes is still unclear. 
Through such bioinformatic tools as PicTar, miRanda, 
and TargetScan, 32,33 target genes of these miRNAs were 
preliminarily screened. The results showed that the target 
genes are those closely related to proliferation, cell-cycle 
regulation, and signal transduction, which are involved in 
aspects of tumorigenesis, tumor development, and metastasis. 
Bioinformatic computational methods for predicting target 
genes are necessary to guide experimental validation. 

One limitation of this study was the lack of cell-function 
experiments and target-gene function validation. Whether 
the abnormally expressed miRNAs identified by microarray 
and verified by qRT-PCR are involved in the pathogenesis of 
laryngeal carcinoma requires further verification by function 
experiments. 

In summary, miRNA-expression profiling of laryn- 
geal carcinoma was performed using microarray analysis, 
and qRT-PCR verification was conducted. The decreased 
expression of let-7f-5p and miR-195-5p is a novel find- 
ing in head and neck cancer. In addition, target genes of 
identified miRNA were predicted using multiple databases. 



The identified miRNAs may have potential value as novel 
biomarkers and treatment targets of laryngeal carcinoma. 
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